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Introduction. 


Planar  Transmission  Media 

Planar  transmission  media  at  millimeter  wavelengths 
provide  a  reasonably  good  performance  and  lend  themselves  to 
mass  production  techniques.  There  are  three  categories  under 
which  such  lines  may  be  placed;  planar  and  quasi-planar  (i.e. 
microstrip  line  and  its  variations,  and  finline),  dielectric  guides 
(dielectric  slab,  image  ‘line,  insular  guide,  inverted  image  line,  etc.), 
and  H-guides  (groove  guide,  trough  guide,  etc.)  [14]. 

At  millimeter  wavelengths,  the  most  commonly  used 
planar  transmission  lines  are  microsirip  and  microstrip-like 
(inverted  and  suspended)  lines.  These  lines  are  suitable  for  the 
design  of  low-cost,  mass-producible  millimeter  wave  integrated 
circuits.  For  microwave  integrated  circuits  up  to  110  GHz,  the 
main  contenders  are  microstrip,  suspended  microstrip,  fin  line  and 
image  line.  Some  examples  of  these  lines  are  presented  in  figure 
1. 


Microstrip  and  microstrip-like  transmission  lines 
consist  mainly  of  a  thin  strip  conductor  on  a  homogeneous  or 
inhomogeneous  dielectric  substrate  that  is  backed  by  a  ground 
plane  of  infinite  conductivity.  Many  numerical  and  analytical 
techniques  exist  which  are  used  to  analyze  the  behavior  of  these 
media  [15].  One  of  the  simplest  of  these  methods  is  the  quasi¬ 
static  approach  [6].  This  approach,  however,  has  a  limited  range 
of  validity,  as  the  nature  of  the  mode  of  propagation  in  this  case  is 
assumed  to  be  pure  TEM,  and  the  transmission  line  characteristics 
are  calculated  from  the  electrostatic  mutual-  and  self  capacitances 
and  inductances  of  the  structure.  The  quasi-static  analysis  is 
therefore  adequate  for  designing  circuits  only  when  the  strip 
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width  and  the  substrate  thickness  are  very  small  compared  to  the 
wavelength  in  the  dielectric  material. 
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The  full-wave  approach,  on  the  other  hand,  is 
complete  and  rigorous.  It  takes  into  account  the  hybrid  nature  of 
the  mode  of  propagation,  and  the  transmission  line  characteristics 
are  calculated  by  determining  the  propagation  constant  of  the 
device.  The  aforementioned  hybrid  modes  are  a  superposition  of 
TMy  and  TE^  fields  that  may,  in  turn,  be  expressed  in  terms  of  two 

scalar  functions,  Vjf  and  \|/,  respectively. 


There  are  several  methods  available  for  calculating  the 
propagation  constant,  including  the  integral  equation  methods, 
finite  difference  methods,  spectral  domain  methods  [2,8,18],  and 
mode  matching  [7,10,11,19]. 
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Figure  1.  Some  examples  of  planar  transmission 
lines;  (a)microstrip  line,  (b)inverted 
microstrip  line,  (c),  (d)coplanar 
waveguide,  and  (e)coupled  microstrip 
line. 


Chapter  1.  Mode  Matching. 

Mode  matching  is  one  of  the  most  frequently  used 
techniques  for  formulating  boundary  value  problems.  Its  primary 
advantage  is  that  it  does  not  generate  spurious  solutions.  It  does, 
however,  have  a  relative  convergence  problem,  so  the  accuracy  of 
the  results  should  be  verified  carefully.  It  is  not  a  very  efficient 
method,  and  not  suitable  for  CAD  packages,  but  it  does  afford  an 
exact  solution  within  a  logical  margin  of  error. 

Mode  matching  is  used  when  the  structure  in  question 
can  be  identified  as  the  junction  of  two  or  more  regions,  each 
belonging  to  a  separable  coordinate  system.  For  a  planar 
structure,  rectangular  cartesian  coordinates,  which  are  a  separable 
system,  are  used  to  describe  the  structure.  One  other 
consideration  in  mode  matching  is  that  in  each  region,  there  must 
exist  a  set  of  well-defined  solutions  of  Maxwell's  equations  which 
satisfy  all  the  boundary  conditions  of  the  structure,  except  the 
continuity  conditions  at  the  junctions  between  the  regions.  Thus, 
the  separation  of  the  structure  into  regions  must  be  done  in  a 
well-defined  and  judicious  manner  such  that  a  simple  and  well¬ 
converging  solution  may  be  obtained. 

The  steps  followed  in  the  mode-matching  procedure 
are  simple  and  straightforward.  The  first  step  is  to  define  a 
certain  set  of  normal  basis  functions  for  each  region  of  the  device, 
and  to  expand  the  unknown  fields  in  these  regions  with  respect  to 
these  normal  functions.  For  the  sake  of  simplicity,  these  basis 
functions  will  be  called  "modes",  although  they  do  not  satisfy  the 
source-free  wave  equation  with  all  the  boundary  conditions.  They 
do,  however,  satisfy  the  wave  equation  in  their  respective  regions, 
and  they  will  be  subject  to  the  boundary  conditions  of  those 
regions.  The  functional  forms  of  these  modes  are  already  known, 
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so  the  electric  fields  are  now  actually  defined  by  the  weight  of 
each  mode.  In  this  manner,  the  original  problem  reduces  to  that 
of  determining  the  set  of  modal  coefficients  associated  with  the 
field  expansions  in  the  various  regions.  This  procedure  leads  to  an 
infinite  set  of  linear  simultaneous  equations  for  the  unknown 
modal  coefficients.  To  obtain  the  exact  solution  to  the  problem  at 
hand,  one  must  solve  this  infinite  set  of  equations,  a  generally 
impossible  task.  Approximation  techniques,  such  as  truncation  or 
iteration,  must  therefore  be  applied,  and  herein  lies  the  difficulty 
of  mode  matching.  In  a  straightforward  analysis,  the  number  of 
modes  retained  will  determine  the  accuracy  of  the  solution.  More 
modes  would  logically  seem  to  provide  a  more  accurate  solution. 
However,  computation  time  increases  as  the  square  of  the  number 
of  modes  retained.  This  is  one  important  factor  which  necessarily 
limits  the  number  of  modes  which  can  be  retained  in  each  region. 
Another  source  of  numerical  difficulties  arises  from  the  fact  that 
planar  structures  such  as  the  microstrip  have  geometrical 
discontinuities  in  the  form  of  sharp  edges.  In  this  case,  the  fields 
must  be  subjected  to  one  more  physical  condition,  known  as  the 
edge  condition,  which  states  that  the  power  of  the  electric  and  the 
magnetic  fields  at  the  edges  must  be  finite.  This  extra  condition  is 
needed  so  that  a  unique  solution  to  Maxwell's  equations  may  be 
obtained.  In  mode  matching,  this  condition  translates  into  a 
relationship  between  the  number  of  modes  retained  in  each 
region  and  the  actual  physical  dimensions  of  the  regions.  This 
relationship  will  become  more  clear  as  the  analysis  of  the 
particular  device  under  consideration,  the  coplanar  microstrip 
line,  progresses. 

The  mode  matching  technique  may  be  extended  to 
include  cases  of  continuous  spectra.  This  work,  however,  will  not 
take  such  a  case  into  consideration. 


Chapter  2.  Analysis  of  the  Coupled  Microstrip 

Line 

2.1.  Parallel-plate  Waveguides 


Figure  2.  Parallel-plate  waveguide.  The  medium 
between  the  plates  is  described  by  the 
general  dielectric  constant  e(x,y). 

A  parallel-plate  waveguide  (figure  2)  is  a  rectangular 
structure,  so  it  can  easily  be  described  using  rectangular 
cartesian  coordinates.  Without  loss  of  generality,  it  can  be 
assumed  that  the  metal  plates  of  the  waveguide  lie  parallel  to  the 
plane  of  two  of  the  axes,  and  that  the  direction  of  propagation  lies 
on  one  of  these  axes.  The  medium  between  the  plates  may  be 
homogeneous  or  inhomogeneous.  In  the  case  of  a  homogeneous 
medium,  finding  the  field  distribution  is  elementary  and  the 
procedure  may  be  found  in  any  textbook.  In  the  latter  case  of  an 
inhomogeneous  medium,  the  case  where  the  medium  is  layered 
in  only  one  direction,  parallel  to  the  plane  of  the  metal  plates  of 
the  waveguide,  will  be  investigated.  Assuming  that  the  plates  lie 
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parallel  t<"  the  x-z  plane,  and  that  the  direction  of  propagation  is 
along  the  z  axis,  TMy  and  TEy  solutions  to  Maxwell’s  equations 
may  be  constructed.  Using  untilde’d  variables  to  denote  TMy 
quantities,  and  tilde'd  variables  to  denote  TEy  quantities,  the  TMy 
field  components  may  be  written  as  (see  Appendix  B)  [13]: 


1  3  f(x,y)  -jk^; 


jo)e(y)  ax  ay 


1 


Ev  = 

'  J  w  e(y) 


ay^ 


f(x,y)  e 


-jk^ 


Hx  =  jkzf(x,y)e 


-jki' 


Hy=0 


(2.1.1) 


^  _  k^  af(x,y)  _-jk^ 
CO  e(y)  ay 


Hz  = 


af(x,y)  _.jk^ 
- e 

ax 


where  f(x,y)  is  the  TMy  scalar  potential,  e(y)  is  the  y-dependent 
permittivity,  k^=(0^|J.o£(y)  Is  the  wavenumber,  and  is  the 

propagation  constant  in  the  z  direction  [4,7,19]. 


For  the  TEy  field  components,  the  corresponding 
equations  are: 


Ex  =  -j  k  J(x,y) 


E 


2 


5f(x,y)  -jk/ 
- e 

dy 


2^ 

1  a  f(x,y)  -jk/ 
- e 

j  to  Po  dx  dy 


Hy  = 


1 


jtOPo 


.2  a 

k  +■ 


2  \ 


dy^l 


f(x.y)  e 


■jk^ 


kj  af(x,y)  -jk^ 
- e 

copo  dy 


(2.1.2) 


Here,  f  (x,y)  is  the  TEy  scalar  potential. 
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In  a  parallel  plate  waveguide  with  the  x  dimension  of 
the  waveguide  very  large,  or  for  high  frequencies,  the  potential 
functions  for  the  TMy  and  the  TEy  fields  may  be  assumed  to  be  a 

function  of  y  only.  So,  the  TMy  potential  may  be  written  as: 

f(x,y)  =  V(y)  (2.1.3) 

and  the  TEy  potential  function  may  be  written  as: 

f(x,y)  =  \i/(y)  (2.1.4) 

Substituting  (2.1.3)  into  equations  (2.1.1)  for  the  TMy 
fields  will  yield: 

Hx  =  jk,V(y)e'''^' 

Hy  =  0  (2.1.5) 

H,=  0 

Similarly,  substituting  (2.1.4)  into  equations  (2.1.2) 
for  the  TEy  field  components  will  yield: 


Ex  =  0 


T-  1  2  .  -ikz 

E.,  = -  k  +  — -  v(y)  e  •' 

2 


j«e(y)\  dy^ 


kj  dY(y).,-jk? 
- ; — e 


CO 


e(y)  dy 


Ex  =  'j  Jf/  V(y)e 


=  0 


Ey=0 
E,=  0 


Hy  = 


1  2  d 
k  + 


dy^/ 


\j/(y)e 


(2.1.6) 


~  d\j/(y)^.jk^ 

"z - 1 - S 

copo  dy 


2.2. 


Three-Layer  Parallel-Plate  Waveguide: 


Figure  2.2.1.  Three-layer  parallel  plate  waveguide. 

At  this  point,  it  would  be  advisable  to  separate  the 
structure  into  regions.  For  a  three-layer  parallel-plate  waveguide 
(a  parallel  plate  waveguide  with  three  layers  of  dielectric),  the 
structure  can  be  separated  into  three  regions,  each  region 
characterized  by  its  own  dielectric.  So,  for  the  parallel-plate 
waveguide  in  fig(2.2.1),  region  1  will  be  that  region  where  0:^<hi, 
h  1  being  the  point  where  the  dielectric  constant  of  the  medium 
changes  from  E]  to  £2-  Region  2  will  be  that  region  where  hi^<h2, 
h2  being  the  point  of  transition  from  £=£2  to  £=£3.  Similarly,  region 
3  will  be  that  region  for  which  h2^^h,  h  being  the  height  where 
the  top  metal  plate  of  the  waveguide  is  situated.  Each  region  will 
be  characterized  by  its  own  set  of  TM^  and  TEy  potential 
functions,  f(y)  being  the  TMy  potential,  and  f  (y)  the  TEy  potential 
of  region  i,  i=l,2,3.  Both  the  TMy  and  TEy  scalar  potentials  are 
solutions  of  the  scalar  Helmholtz  equation: 
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— x(y)  =  0 

dy  (2.2.1) 

a  general  and  complete  solution  to  the  above  equation  being; 

X(y)  =  tti  cos  (Ky)  +  a2 sinf  (Ky)  (2.2.2) 

where  aj  and  02  are  constants  which  depend  on  the  boundary  and 
continuity  conditions  of  the  structure. 

The  function  sinf(Ky)  is  defined  as  sinf(Ky)=sin(Ky)/K 
and  its  usefulness  will  be  presented  later  on  in  this  analysis. 

The  electric  and  magnetic  fields  must  conform  to 
certain  boundary  and  continuity  conditions.  These  conditions  will 
serve  to  specify  the  unknown  variables  ai,a2.and  Kin  equation 
(2.2.2).  The  general  solution  to  the  scalar  wave  equation  shown 
above  is  modified  such  that  these  continuity  conditions  can  be 
implemented. 

Substituting  x(y)  into  the  TM^  equations,  and 
applying  the  conditions  of  zero  tangential  electric  field  on  the  pec 
plates  will  yield  the  general  form  for  the  potential  for  the  TM^ 
modes  of  a  three  layer  parallel'plate  waveguide; 


a  cos  ( kyi  y ) 


;0.<y  <hi 


¥(y)  = 


b  cos[  ky2(h2-y)] 

+  c  sinf[  ky2(h2-y)] 

d  cos[  ky3(h-y)] 


;hi<y<h2  (2.2.3) 


;h2Sy<h 


Correspondingly,  substituting  the  appropriately 
modified  forms  for  the  corresponding  form  for  the  TEy  scalar 
potential  into  the  TEy  field  equations  will  yield; 


a  sinf  ( kyi  y ) 


;0.<y<hi 


b  cos[  ky2(h2-y)] 

+  c  sinfl  ky2(h2-y)] 

d  cos[  ky3(h-y)] 


;hi<y<h2  (2.2.4) 


;h2<y<h 


In  general,  the  coefficients  of  the  TMy  and  TEy  scalar 
potentials,  as  well  as  their  corresponding  eigenvalues,  will  not  be 
identically  equal. 

The  continuity  conditions  between  the  different 
regions  state  that  the  tangential  fields  at  a  dielectric  discontinuity 
must  be  continuous.  Applying  these  conditions,  and  solving  the 
resulting  equations,  will  yield  the  eigenvalue  equations  for  the 
TMy  and  TEy  y-directed  eigenvalues  as  well  as  the  corresponding 
expansion  coefficients  of  the  scalar  potentials  in  equations  (2.2.3) 
and  (2.2.4).  The  system  is  underdetermined,  so  the  four 
expansion  coefficients  in  each  case  will  be  found  within  a 
multiplicative  constant. 

Following  the  procedure  described  above  will  yield 
the  TMy  eigenvalue  equation  for  the  three-layer  parallel-plate 
waveguide; 


(2.2.5) 


^^tan(kyihi)  +  ^^tan[ky2(h2-hi)]  +  ^tan[ky3(h-h2)] 
Ej  62  £3 

tan(ky ih  1)  tan[k  2(h2-h  1)]  ^  tan[ky,(h-h2)] 

Ej 

and  the  coefficients  are: 


1 .  ky3  real 

n.  /i,  V,  M  ^aginary 

cos[ky3(h-h2)] 

(2.2.6. a) 

E-  .  . 

-  —  ky3  sin[ky3(h-h2)]  ky3  real 

£3 

-  —  ky3  tan[ky3(h-h2)]  ky3  imaginary 

£3 

(2.2. 6.b) 

cos[ky3(h-h2)]  ky3  real 

1 .  ky3  imaginary 

(2.2.6.C) 

b  cos[ky2(h2-hi)]  +  c  sinlTky2(h2-h  j)] 
cos(kyihi) 

(2.2.6.d) 

One  other  relationship  of  great  importance  in  this 
analysis  is  the  dispersion  equation.  This  equation  links  the 
eigenvalues  k^,  ky,  and  k^,  with  the  wavenumber  k  in  a  certain 
medium.  The  dispersion  equation  states  that: 

2222  2 

k,  +  ky  +  kz  =  k  =0)  P(^(y)  (2.2.7) 


for  the  TMy  case,  and  correspondingly: 


(2.2.8) 


1  2 


~2  ~2  2  2  2 
kx  +  ky+k2  =  k  =0)  fi(^(y) 


for  the  TEy  case. 

For  the  particular  case  of  the  three-layer  parallel- 
plate  waveguide,  it  has  been  assumed  that  kxi=0,  i  =  1,  2,  and  3, 

so  the  dispersion  relations  for  the  TMy  and  the  TEy  case 
correspondingly  become: 

2  2  2  2 

kyi  +  kz  =  kj  =  CD  p,{]pj  (2.2.9) 


and 


~2  2  2  2 

kyi -f- kz  =  kj  =  (0  (2.2.10) 

kz  being,  of  course,  the  propagation  constant  of  the  waveguide, 
and  k  the  wavenumber.  The  z  axis  (and  therefore,  the  direction 
of  the  propagation  constant)  are  parallel  to  the  planes  of  the 
discontinuities,  and  so  continuity  forces  kz  to  be  the  same  in 
every  region.  This  gives  rise  to  a  very  useful  relationship 
between  the  eigenvalues  in  the  y  direction  in  each  region.  This 
relationship  is  derived  by  subtracting  the  dispersion  relation 
defined  in  the  one  region  from  the  corresponding  equation  in  the 
other  region.  Thus,  the  dispersion  relationships  corresponding  to 
two  neighboring  regions  of  the  waveguide  reduce  to: 

2  2  2 

kyi  -  Ho(£i  -  £j)  (2.2.11) 

and: 

~2  ^2  2 
kyi-i-kyj=co  iio(£i-ej) 


(2.2.12) 


so  that  the  y  eigenvalues  in  each  region  of  the  waveguide  are  not 
independent  variables,  but  are  linked  through  the  dispersion 
relation. 


The  TEy  eigenvalue  equation  as  well  as  the 
coefficients  for  the  TEy  scalar  potential  may  be  derived  in  a 
manner  similar  to  that  used  in  the  derivation  of  the  TMy 
potential.  The  TEy  eigenvalue  equation  is  thus  found  to  be; 


tan(kyihi)  ^  tan[ky2(h2-hi)]  ^  tan[ky3(h-h2)] 

kyl  ky2  ky3 

tan(kyihi)  ~  rr  /i.  m  tan[ky3(h-h2)] 
- ky2tan[ky2(hrhi)] — ^ 

kyl  ky3 


(2.2.13) 


and  the  coefficients  of  the  TEy  scalar  potential  in  the  waveguide 
are; 


c  = 


b  = 


1 .  ky3  real 

1  — 

- : - ky3  imaginary 

cos[ky3(h-h2)] 

cos[ky3(h-h2)]  ky3  real 

1 .  ky3  imaginary 

sinf[ky3(h-h2)]  kysreal 

sinflky3(h-h2)]  ~  . 

- ^ ^  ky3imagmary 

cos[ky3(h-h2)] 


(2.2. 14. a) 


(2.2. U.b) 


(2.2.14.C) 


(2.2.14.d) 


~  _  b  cos[ky2(h2-hi)]  +  c  sinf[ky2(h2-h  j)] 

3. - zz - 

cos(  kyih^ ) 


Working  in  the  reverse  order  now,  the  eigenvalue 
equation,  in  conjuction  with  the  modified  dispersion  relation,  will 
yield  the  eigenvalues  in  the  y-direction  of  the  three-layer 
parallel-plate  waveguide.  The  dispersion  relation  will  then  yield 
kz,  the  propagation  constant  of  the  waveguide  in  the  z  direction: 


/  2  ~2 
k2=V(0  |idEi-kyi 

for  the  TMy  case,  and 

k2=  y  (0  ^idEi-kyi 
for  the  TEy  case. 


(2.2.15) 


(2.2.16) 


In  addition,  the  ky,  and  kyj  eigenvalues  serve  to 
specify,  within  a  multiplicative  constant,  the  expansion 
coefficients  of  the  corresponding  scalar  wave  equations.  These 
scalar  wave  equations,  substituted  into  the  corresponding  TMy  or 
TEy  field  equations  will  yield  the  TMy  or  TEy  electric  and 
magnetic  field  components.  Since  the  coefficients  are  found  only 
within  a  multiplicative  constant,  only  the  distribution  of  the  E 
and  the  H  fields  in  the  waveguide  can  be  found. 

When  the  TMy  and  TEy  scalar  potentials  were  defined, 
a  new  functional  form  was  introduced.  This  form  was  the 
function  sinf(Kx),  which  was  subsequently  defined  as: 

sinf(Kx)  = 

K  (2.2.17) 
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From  the  definition  of  the  "sinf"  function  as  stated 
above,  it  is  evident  that  this  function  has  the  following 
properties: 

^  {sinf(Kx)}=cos(Kx) 

ay  (2.2.18) 

and 


lim 

ky  — >0 


{sinf(Kx))=  X 


(2.2.19) 


The  "sinf"  function  is  preferred  over  the  simpler  "sin"  because  it 
offers  two  very  important  advantages. 

One  advantage  is  that  sinf(Kx)  offers  the  correct 
solution  to  the  scalar  Helmholtz  equation  for  vanishing  K.  Using 
sin,  the  general  solution  to  the  scalar  wave  equation  is: 

X(y)  =  a  1  cos(Ky)  +  ai  sin(Ky)  ^  2  201 


As  K  approaches  0,  x(y)  becomes: 
X(y)  =  (Xi 


(2.2.21) 


Although  this  does  provide  a  solution,  it  is  not  a  general  solution 
to  Laplace's  equation. 


With  the  "sinf"  representation,  however,  the  solution 
to  Laplace's  equation  becomes: 


X(y)  =  oci  cos(Ky)  +  a2  sinf(Ky) 


(2.2.22) 


which,  for  vanishing  K  yields: 

X(y)  =  ai  +a2y 


(2.2.23) 


This  latter  form  does  provide  a  general  solution  to  Laplace's 
equation. 


The  second  very  important  advantage  of  the  "sinf" 
function  over  the  "sin"  is  found  in  the  nature  of  the  eigenvalues 
K.  For  a  purely  lossless  dielectric  medium,  K  will  be  either  real  or 
imaginary.  For  K  imaginary,  the  "sin"  function  would  also  be 
imaginary,  which  fact  might  give  rise  to  a  complex  analysis. 
However,  the  "sinf"  function  remains  real  irrespective  of  whether 
K  is  real  or  imaginary.  This  provides  a  distinct  improvement 
over  the  "sin". 


Here,  the  asymptotic  expression  for  "sinf"  is  used  for  the  case 
iKxkO.l. 


2.3.  Four-layer  parallel-plate  waveguide 


Figure  2.3.1.  Four  layer  parallel-plate  waveguide. 

The  TMy  scalar  potential  for  the  case  of  the  four-layer 
parallel-plate  waveguide  may  be  written  as; 


V(y)  = 


a  cos(kyiy) 

b  cos[ky2(hry)] 

•f  c  sinf[ky2(h2-y)] 
d  cos[ky3(h3-y)] 

+  e  siiif[ky3(h3-y)] 


;  0  <  y  <hi 


;  hi<  y  <h2 


;h2<y  <h3 


(2.3.1) 


f  cos[ky4(h-y)] 


;  h3  <  y  <  h 


The  eigenvalue  equation  is  found  using  the  same  procedure 
that  was  described  for  the  three-layer  parallel-plate  waveguide. 
The  TMy  eigenvalue  equation  is  therefore: 

^  tan(kyihi)  +  ^  tan[ky2(h2-hi)] 

El  £2 

-I-  ^  tan[ky3(h3-h2)]  +  ^  tan[ky4(h-h3)]= 

£3  £4 

^  tari(kyihi)  tan[ky2(h2-hi)]  ^  tan[ky3(h3-h2)]  + 

£1  ky2  £3 

^  tan(kyihi)  ^  tan[ky2(hrhi)]  ^  tan[ky4(h-b3)]  -h 
£1  "  ky2  -  £4 

^  tan(kyihi)  tan[ky3(h3-h2)]  ^  tan[ky4(h-h3)]  -h 

£1  ky3  '  £4 

tan[k.:(h2-h,)l  p-  tan[kj,3(h3-h2)l  ^  tan[ky4(h-h3)] 

£2  k,3  £4  ,2.3.2) 

The  coefficients  which  serve  to  define  the  TM^  scalar  potential  of 
the  waveguide  are,  within  a  multiplicative  constant: 


1 .  ;  ky4  real 

n,  ^  u  M  ’ 
cos[ky4(h-h3)] 

(2.3.3. a) 

£, 

ky4  sin[ky4(h-h3)]  ;  ky4  real 
£2 

£, 

■—  ky4  tan[ky4(h-h3)]  ;  ky4  imaginary 
£2 

(2.3.3.b) 

cos[ky4(h-h3)]  ;  ky4  real 

1.  ;  ky4  imaginary 

(2.3.3.C) 
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tan(kyihi)  taii[ky2(h2:-hi)] 

- == -  T - = - 

kyl  ky2 

^  taii[ky3(h3-h2)]  ^  tan[ky4(h-h3)]  _ 
ky3  ky4 

(2.3.5) 

k,,ian[k,^h^h,)l  ‘^tMh-hs)l  ^ 

kyl  ky4 

!g?gy'!l->.^3an[Mhrh.)l  ‘^iMh-h3)l  ^ 

kyl  ky4 

i?M2Mk,3ttn[Mhrh,)l 

ky2  ky4 

The  coefficients  of  the  scalar  potential  are,  within  a  multiplicative 
constant: 


1. 

;  ky4  real 

1. 

;  ky4  imaginary 

cos[ky4(h-h3)] 

(2.3.6. a) 

cos[ky4(h-h3)] 

;  ky4  real 

1. 

;  ky4  imaginary 

(2.3. 6.b) 

sinf[ky4(h-h3)] 

;  ky4  real 

tan[ky4(h-h3)] 

;  ky4  imaginary 

(2.3.6.C) 

b  =  d  cos[ky3(h3-h2)]  +  e  sinflky3fb3-h2)] 


(2.3.6.d) 
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c  =  -d  ky3  sin[ky3(h3-h2)]  +  e  cos[ky3(h3-h2)] 

~  _  b  cos[ky2(h2-hi)]  +  c  sinf[ky2(h2-hi)] 

3. - :::::: - 

cos(kyihi) 

Thus,  so  far,  the  TM^  and  TEy  potential  functions  for 
the  three-  and  four-layer  parallel-plate  waveguides  have  been 
derived.  Their  connection  with  the  coupled  microstrip  line,  which 
is  the  main  subject  of  this  paper,  will  become  apparent  shortly. 


(2.3.6.e) 

(2.3.6.f) 


Chapter  3.  Coupled  Microstrip  Lines 

Coupled  microstrip  lines  are  used  in  a  number  of 
circuit  applications,  principally  as  directional  couplers,  filters,  and 
delay  lines.  Mode  matching  will  be  used  to  calculate  the 
propagation  constant  of  such  a  line. 

A  pair  of  microstrip-like  transmission  lines  as  shown 
in  fig.  (3.1)  are  known  to  have  the  property  of  a  broadband 

directional  coupler  when  placed  in  parallel  proximity  to  each 
other.  As  a  result  of  this  proximity,  a  fraction  of  the  power 

present  on  the  main  line  is  coupled  to  the  secondary  line.  The 

power  coupled  is  a  function  of  the  physical  dimensions  of  the 
structure  and  the  direction  of  propagation  of  the  primary  power. 


Figure  3.1.  The  coupled  microstrip  line. 

In  general,  the  coupled  line  structures  shown  in  fig. 
(3.1)  support  two  modes,  the  even  and  the  odd  modes.  This 
induces  the  coupling  between  the  two  transmission  lines.  The 
properties  of  the  coupled  structures  may  be  described  in  terms  of 
a  suitable  linear  combination  of  these  even  and  odd  modes.  For 
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geometrically  symmetric  structures,  the  geometrical  plane  of 
symmetry  may  also  be  thought  of  as  an  electrical  plane  of 
symmetry.  The  even  mode  has  equal  amplitude  and  equal  phase 
correspondence  with  respect  to  this  plane  of  symmetry,  in  which 
case  this  plane  takes  the  form  of  an  open  circuit.  In  the  even 
mode  case,  therefore,  the  plane  of  symmetry  takes  the  form  of  a 
perfect  magnetic  wall  (pmc).  The  odd  mode,  on  the  other  hand, 
has  an  equal  amplitude  180®-out-of-phase  correspondence  with 
respect  to  the  plane  of  symmetry,  so  this  plane  acts  like  a  short 
circuit  and  is  simulated  by  a  perfect  electric  conductor  (pec)  wall. 
Thus,  the  directional  coupler  can  be  treated  as  a  two-port 
network,  and  the  total  response  can  be  obtained  by  superimposing 
the  responses  calculated  for  the  even  and  odd  mode  excitations. 
This  reduces  the  problem  to  about  half  its  original  size,  which  is  a 
great  advantage  of  symmetric  structures. 


3.1.  Propagation  Constant 

The  structure  shown  in  figure  (3.1)  is  symmetric  with 
respect  to  a  plane  drawn  perpendicular  to  the  x-y  and  x-z  planes 
and  placed  halfway  between  the  two  metal  strips.  Thus,  the 
structure  under  consideration  is  modified  to  that  shown  in  figure 
(3.1.1). 


Figure  3.1.1.  Symmetry  as  applied  to  the 
coupled  microstrip  line 

As  was  mentioned  earlier,  the  left-most  boundary 
becomes  either  a  pmc  or  a  pec  wall,  according  to  whether  the 
search  is  conducted  for  the  odd-mode  or  even-mode  propagation 
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constant.  There  is  no  perpendicular  boundary  to  the  right  of  the 
microstrip,  so  propagation  in  the  positive  x  direction  for  large  x 
must  take  the  form  exp(-jkxx),  such  that  there  is  an  exponential 

decay  for  x— This  serves  to  put  a  lower  bound  on  the  search 
for  k2,  since  this  exponential  decay  is  observed  only  for  the  case 
where  k^  is  imaginary.  The  dispersion  relation  states  that: 

2  2  2  2  2 

kx+ky+k2  =  k  =  0)  (3.1.1) 

so,  for  kx: 

2  2  ,22 

kx  =  CO  PqE  -ky-k^  (3.1.2) 

The  condition  for  imaginary  kx  yields:  _ 

(0  Poe-ky-kz<0  (3.1.3) 


or: 


kz> 


V 


2 

to  PdE- 


(3.1.4) 


An  upper  bound  may  also  be  placed  on  k^,  since  it  may 
never  be  larger  than  the  unbounded  propagation  constant  in  the 
medium  with  the  highest  relative  dielectric  constant.  Thus: 

k^<  (oVp(^(^ax(ei)  (3.1.5) 


The  structure  must  further  be  separated  into  regions 
for  modal  expansion.  Two  more  bounding  planes  may  be  drawn 
parallel  to  the  plane  of  symmetry,  at  each  end  of  the  metal  strip. 
Also,  since  the  structure  must  be  bounded,  a  top  plate  must  be 
placed  above  the  structure  at  such  a  distance  that  it  does  not 
interfere  with  the  computation  of  the  propagation  constant  of  the 
actual  open  structure  (this  top  plate  was  also  included  in  figure 
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3.1.1).  The  two  lateral  bounding  planes  may  be  thought  of  as 
perfect  magnetic  conductor  walls,  since  the  tangential  electric  and 
magnetic  fields  must  be  continuous  across  them.  The  original 
structure  has  thus  been  separated  into  four  distinct  regions, 
shown  in  figure  (3.1.2). 


Figure  3.1.2.  The  four  distinct  regions  of  the 
coupled  microstrip  line 

As  was  stated  earlier,  a  set  of  normal  modes  must  be 
defined  for  each  region,  and  the  unknown  fields  in  each  must  be 
expanded  with  respect  to  these  modes.  Each  region  may  be 
thought  of  as  a  multilayered  parallel-plate  waveguide,  consisting 
of  two  parallel  metal  sheets  placed  at  a  certain  distance  from,  and 
parallel  to,  each  other,  with  a  layered  dielectric  medium  between 
them.  Thus,  regions  1  and  4  are  four-layered  parallel-plate 
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waveguides,  region  2  is  a  three-layered  parallel -plate  waveguide, 
and  region  3  is  a  one-layered  parallel-plate  waveguide.  This 
latter  case  has  been  studied  extensively  in  the  literature,  and  has 
not  been  presented  here.  The  normal  modes  for  the  three-  and 
four-layered  parallel-plate  waveguides  have  been  defined  in 
chapters  2.2  and  2.3,  respectively.  It  now  remains  to  establish  the 
equations  which  define  the  electric  field  components  with  respect 
to  these  modes.  Let  it  be  noted  that  the  following  modal 
expansions  and  field  representations  are  not  unique,  but  just 
constitute  one  possible  solution. 

Using  the  expressions  derived  for  the  TMy  and  TEy 
modes  of  the  three-  and  four-layer  parallel-plate  waveguides,  the 
corresponding  expressions  for  the  electric  and  magnetic  field 
expansions  may  be  derived.  As  shown  in  Appendix  B,  the  electric 
and  magnetic  fields  in  each  region  due  to  the  TMy  mode  potentials 
are  given  by: 


E  •  =  — 

^X1  ^ 


3  (fi(x,y)e'-’*^^ 

d\  dy 


Hxi  = 


3(/i(x,y)e’-’’'^ 


E  ■  =  ■ 


yi  ax^ 


+  kf  (fi(x,y)e'-'‘'^ 


^  1  3  (Ji(x,y)e'-'‘'^ 

-V 


dy  dz 


Hyi=0 


H^i  = 


(Ji(x,y)e''‘'^ 


(3.1.6) 


and  the  electric  and  magnetic  fields  due  to  the  TMy  fields  are: 
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Exi  = 


_  3(fi(x.y)e'-’^^ 
dz 


Eyi  =  0 


Ezi=- 


a(Ji(x.y)e'^‘^^ 
3x 


Hxi  = 


1  3  (/i(x,y)e‘-’’‘^) 


z 

/  2 


3x  9y 


Hyi=^ 


+  k 


ax' 


(j’i(x,y)e'-’’‘^  (3.1.7) 


Hzi  = 


1  a\ji(x,y)e'-'’‘^ 
z  ay  az 


where  yi  =  jcoe,  and  z  =  jco^oll^],  and  subscript  i=l,2,3,or  4, 
correspondingly,  for  each  region. 


The  electric  and  magnetic  potentials  in  each  region  are; 


Region  1: 

f  im(x,y)  = 

(3.1.8.a) 

Region  2: 

f2m(x,y)  =  V2m(y)^  A2mCOS[k^2m(x-Wi)] 

+  B2mSinf[k^2m(x-Wi)l} 

/2m(x,y)  =  V2m(y)  {  A2mCOS[k,2m(x-W  i)] 

- .  +  B  2n,sinf[k^2m(x-w  i)]  } 

Region  3: 

/3m(x,y)  =  V3m(y)  (  A3nfOS[k;^3 Jx-W,)] 

+  B3^inf[kx3m(x-Wi)]} 
f 3m(x,y)  =  V3m(y)  {  A3mC0S[k,3jx-Wi)] 

+  B3mSinf[k^3jx-Wi)]} 

Region  4: 

J4m(x,y)  =  V4in(y)  exp[-jk4^m(x-W2)]  (3. 1 .8.d) 

/4m(x,y)  =  V4m(y)  exp[-jk4^  Jx-Wj)] 

where  \)/i(y)  and  Vi(y)  are  the  TM^  and  TE^  potential  functions  for 
the  corresponding  cases  of  the  one-,  three-,  and  four-layer 
parallel-plate  waveguides. 

For  a  suppressed  z-dependence,  the  equations  for  the 
electric  and  magnetic  field  components  yield,  for  each  region: 


(3.1.8.b) 


(3.1.8.C) 


REGION  1 

Electric  Field  Components: 


1 


jcoei(y)  *"=1 


dy 


-k 


Ixm 


Sin 

cos 


(klxmX)Ai 


m 


M,  ^ 

-jkzS  Vim(y) 

m=l 


M, 

jOEi(y)  m=l 


1 


Iz- 


jc«i(y) 

M,  ^ 


.  v-.dViiy)cos 

ivyi  m=l 


m 


(3.1.9.a) 


Magnetic  Field  Components: 

Ml 


m=l 


1  v'  ^Vim(y) 


jOjPo  m=l 


dy 


1  Icos.r  \  7 

~  (^Ixnn?^)  ^Im 

\-klxm/^‘" 


=  -^'L  (kz+kixm)¥lm(y)^^"f(klxmX) 


ly 


”>=1 


sinf^r 
cos 


V*  /  'il 

Hlz=Z¥lm(y)  , 

m=l  \  ^  I 


Sin 

cos 


(kixinX)Ai, 


0)pQin=l 


dy  cos 


(3.L9.b) 
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REGION  2 

Electric  Field  Components: 

jcoe2(y)  m=i 

+  B2niCOS[(k2xm(x-W,)]}  + 

M2  _  _ 

jkzX  V2m(y)  {A2mCOS[k2xm(x-Wl)] 

m=l 

+  B2mSinfl(k2xm(x-Wi)]) 


IVI2 

E2y  = -  X  (kz+kLm)V2n,(y)  t  A2;„COS[(k2xn,(x-Wi)l 

jo3e2(y)  m=l 

+  B2mSinf[(k2xn,(x-w,)3) 

-ky  d\|r-_(y) 

E2z= - X  ~-aZ—  (  A2mCOS[(k2xm(x-Wi)] 

Ci)e2(y)  m=i 

+  B2n^inf[(k2xm(x*Wi)]  }  - 

X  V2m(y)  (  -A2mk2xn,sin[(k2xm(x-Wi)] 

in=l 

+  B2nfOS[(k2xm(x-Wi)]} 


(S.l.lO.a) 
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Magnetic  Field  Components: 

M, 

^2x  =  K  S  Vlm(y)  {  A2„pOS[k2xm(x-Wi)] 
m=l 

+  B2n^infrk2xm(X'W,)]  }  + 

M  ^ 

1  V  ^V2m(y)  ,  7  r  r~  /•  M 

-  L  — - {  *A2mk2xrrPOS[k2xm(x-Wi)] 

30^0  m=l  ^y 

+  B2n^inf[k2*m(x-wi)]  ) 


(S.l.lO.b) 


1  ^  2  ~  ~  ~ 

^2y=-. -  Z  (kz+k2xtn)¥2m(y){  A2mCOS[k2xm(x-W  i)] 

J<^0  '*1=1 

+  B2n,sinf[k2zm(x-wi)]  ) 

Mj 

H2z=  Z  V2m(y)  (  -A2n^2xmSin[k2xm(x-Wi)] 
m=l 

+  B2n,COS[k2;^Jx-Wi)]  )  - 
k^  dV2m(y)  r  7  r  rr  /  X-I 

- L  - aZ - {  -A2n^2xniCOS[k2xm(x-Wi)] 

COHom=l  “y 

+  B2mSinf[k2xn.(x-Wi)]  } 

The  equations  for  region  3  are  identical  to  those  for 
region  2  with  subscript  "2"  substituted  by  subscript  "3". 

REGION  4 


Electric  field  components: 


E4x  =  ^ - ^—r  X  jE4xniexp[-jk4^Jx-W2)]  A4n,  - 

M4  _  _  _ 

jkzX  V4m(y)  exp[-jk4xm(X'W2)]  A4m 
m=l 


3 


M4 


E4y  = 


2  2 


jcoe4(y) 


Y,  (kz+k4xm)V4m(y)  exp[-jk4xn,(x-W2)]  A4m 


O.l.ll.a) 


M. 


^42 - 


dy4m(y) 


co£4(y)  dy 


exp[-jk4 

xm  (X-W2)]  A4n,+ 


M4  _  _  _  _ 

X  ¥4m(y)  jk4xiriexp[(-jk4xm(x-W2)]A4„ 
m=l 


Magnetic  field  components: 


Ml 

H4x  =  jkz  X  '1^4m(y)  exp[k4^„(x-W2)]A4m  - 
m=l 

_J_  Y  jk4xrnexp[-jk4xm(x-W2)]A4f„ 

imo  iS  dy 


H4y  = 


1 


M4 


iTl 


^  (kz+k4xm)V4m(y)exp[-jk4xm(x-W2)]A4m 


O.l.ll.b) 


M4 

H4z=  X  V4m(y)  jk4xmexp[-jk4;,ni(x-W2)]A4ni- 
m=l 


diy^mCy)  exp[-jk4xni(x-W2)lA4in 
“tioiS  <'y 


In  the  above  equations,  the  sums  have  been  truncated 
at  Ml  modes  for  the  fields  in  region  1,  M2  modes  for  those  in 
region  2,  etc.  For  the  solution  of  this  system  of  equations  to  be 


unique,  these  equations  must  form  a  square  matrix,  which  means 
that  : 

2  ^3  j  12) 

The  relationship  in  eq.  (3.1.12)  ensures  that  the 
system  of  equations  will  have  a  unique  solution.  It  does  not, 
however,  ensure  convergence  of  the  system.  This  is  because  one 
very  important  geometrical  factor  has  not  been  considered  yet. 
All  planar  structures  have  discontinuities  in  the  form  of  sharp 
edges.  It  was  noted  earlier  in  this  chapter  that  these 
discontinuities  must  conform  to  the  edge  condition,  that  is,  the 
power  of  the  electric  and  magnetic  fields  at  these  points  must  be 
finite.  This  condition  was  found  to  be  adequately  met  when  the 
number  of  modes  retained  in  each  region  were  such  that: 

Ml _  M4  _  h 

M3  M3  h-h3  (3.1.13) 

and 

Ml  ^  M4  _  h 

M3  M3  h-h3  (3.1.14) 

Thus,  the  edge  condition  will  yield  a  relationship 
between  the  number  of  modes  which  should  be  retained  in  each 
region,  and  their  relative  dimensions. 

Matching  the  tangential  fields  at  the  interface  x=wi 
and  0<y<h3  yields; 
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Eiy(x=Wi)  =  E2y(x=wi)  <=> 

Ml  2  2 

-J—Y  (kz+kixin)¥lm(y),^°V^lxniWi)Aim 

Jt^ei(y),S  ^ 

Ml  2  2 

=  I - 7~~r^.  (kz+kixm)Vlni(y)A2m 

jco£2(y)iri 


Ei2(x=Wi)  =  E2z(x=Wi)  «• 


k,  dVlm(y)cOS,,  K  * 


m 


Ml  ^  /  I 

-X'KimCy 


m=l 


Ixm 


COS 

sin 


(kixmW  i)Ainr- 


M 


M- 


S%^A2„-2v2m(y)B2,n 

^1=1 


®£2(y),S  dy 


Hiy(x=wi)  =  H2y(x=wi)  <=> 

1  Ml  2  ~ 

^ - y  (kz+kixm)¥lm(y)^^^  (kixmWi)Aim 

JW^OiS  ''''' 

^2  2  2 

=  T-^y  (k2+k2xm)y2m(y)A2m 


(3. 1.15. a) 


(3.1.15.b) 


(3.1.15.C) 
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Hi2(x=wi)  =  H22(X=Wi)  <=> 

m=l  V  J 


Ml 


-dW^sind  j 

dv  COS 


dy 


M2 

X  ¥2m(y)B2m 

m=l 


d\l/2m(y)X_ 

- - - ^2m 

dy 


(3.1.15.d) 
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Matching  tangential  electric  and  magnetic  fields  at 
x=w2  for  0<y<h3  yields: 

E2y(x=W2)  =  E4j<X=W2)  <=> 

1  ^..2.1 


X  (k2+k2xm)¥2m(y){A2mCOS[k2xm(W2-Wi)] 


+B2mSinf[k2xm(W2-Wi)]  }  = 

M4 


1 


j“£4(y)S 


_  2  2 

/  .  (kz'^k4xni)M^4in(y)A4 


m 


(3.1.16.a) 


E22(x=W2)  =  E42(x=W2)  <=»  . 

kz  ^  {A2mCOS[k2xm(WrWi)] 


C0E2(y)"i 

+B2n,sinf[k2xn,(WrWi)]  }- 

X  ¥2m(y){-A2mk2xinSin[k2xm(W2-Wi)] 


m=l 


+B2nlCOS[k2xm(W2-W])]  )  = 


- - - ^A4ni+  ^  V4in(y)jk4xmA4m 


coe4(y),S  ^y 


H=1 


(3.1.16.b) 


H2y(x=W2)  =  H4y(x=:W2)  <=> 

2  ^2  2  ~2  ~ 

^  (k2+k2xm) V2m(y)  { A2mtX)S[k2xm(W2-W  1)] 


+B2mSinf[k2xm(WrW,)]  )  = 


M4 


I  •'  2  ~2  ~  ~ 

■; - (kz+k4xin)V4ni(y)A4iTi 


(3.1.16.C) 
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H22(x=W2)  =  H42(x=W2)  «> 

M2 

X  ^2miy)  {-A2mk2xinSm[k2xm(WrWi)] 
ni=l 


+B2mCOS[k2xm(WrWi)]  }- 
-  X  A2rtCOS[k2xin(W2-W,)] 

“^^0  dy 


+B2mSinf[k2xm(W2-Wi)]  }  = 


M4 

X  V4m(y)j^4xmA4m 
m=l 


(O)lo 


dy4m(y) 

^  dy 


n=i 


(3.1.16.d) 


The  equations  for  x=wi  and  x=w2  with  h3<ySh  are 
identical  to  the  corresponding  equations  for  0<y<h3  with  the 
subscript  "3"  substituted  for  subscript  "2". 


As  stated  earlier,  the  orthogonality  of  the  modal 
expansion  functions  Viin(y)  Viin(y)  take  the  form: 


I 


1 


over  one  region ejCy) 


¥im(y)Vm(y)dy  for  m;tn 


(3.1.17) 


and: 


I 


.  ¥im(y)¥in(y)dy  for  m;tn 

over  one  region  (3.1.1o) 


So,  orthogonality  can  be  utilized  here  to  reduce  the 
complexity  of  the  continuity  equations.  By  multiplying  each 
equation  by  the  appropriate  orthogonal  function,  and  integrating 
over  the  appropriate  region  (from  0  to  h3),  the  equations  for 
0<y<h3  may  be  rewritten  as: 


40 


Ml 


2  ,  2  ,,_e2  cos 


m=l 

=  (k^kun)I?nnA2„ 


(3,1 


,  Ml 

_h2  COS  .A 

m=l 


“Z 

CO  “ 


V  t“  / 

+  >  lOrv^rJ 


1 


COS 


(klxmWi)Aini= 


xm 


k,”2  .h2 


h2 


^Ir3W^2m-l4;nB2n 


(0 


m=l 


(3 


/I  ^  NT*'^  sinf  .r  sT 

^  (kz+kixm)l2mn  -  (kixmWi)Aini 


m=l 


COS 

2-2  h2  ~ 

~  (k2+k2xn)l4nnA2n 


m=I  V.  / 


Im 


k,  ^e2  sind 


X  ^3mn  (klxinWi)Aimn 


e2 


Mz 


e2  ~ 


l2nnB2m  T^,^4mnA2m 

®140,S 


(3.1 


(3.1, 


.18.a) 


.l.lS.b) 


.18.C) 


18.d) 


4  1 


(kz+kLn)lSn{  A2rf^OS[k2x„(W2-W  j)] 

+B2nSmf[k2xn(WrWi)]  )  = 

M4 


(3.1.18.e) 


2  2  e2 

/ ,  (kz+k4xni)1 1010^4 


tn 


m=l 


k  h2 

—  y  hmn  {A2mCOS[k2xm(W2-Wi)]+B2mSinf[k2xm(W2rWl)]  ) 

m=l 

l4nn{  -A2iS2xnSm[k2xn(W2-W  i)]+B2nCOS[k2xn(W2-W  j)]  }  = 

M4  h2  h2  ~  ~ 

y  IlmnA4m‘^  y  l2mnjk4xmA4ni 


(3.1.18.f) 


2  —2  h2  ~  ~  ~  ~ 

(kz+k2xn)l4nn{A2nCOS[k2xn(W2-Wi)]+B2nSinf[k2xn(W2-Wi)]  }  = 

2  h2  ~ 

y  (k2+k4xm)l2mnA4m 
m=l 


(3.1.18.g) 


e2 


l2nn  {-A2nk2xnSin[k2xn(W2-Wi)]+B2riCOS[k2xn(W2-Wi)]  } 

e2  ~  ~  ~  - 

- —  y  lLn{A2niCOS[k2xm(WrWi)]+B2mSinf[k2xm(W2-Wi)]  }  = 


(3.1.18.h) 


M4 


e2 


k.  ~ 


y  Ilninjk4xiiiA4in  ■*"  y  ^3mnA4  m 


where: 
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Ilinn=f  '-^viin(y)Vin(y)ciy  =  [  V4m(y)Vin(y)dy 
Jo  e(y)  •'0 

^2mn=\  -7-Vim(y)Vin(y)dy 
•'0  e(y) 

.ci  f”’  1  <i'i'im(y)  ,  ,,  f"’  1  dwy),..  , 

l3mn= - : - Vin(y)‘iy  =  "TT - Z - 

Jo  e(y)  dy  Jo  E(y)  dy 

14^=1  — VLn(y)dy 

Jo  e(y)  d  y 


¥in(y)dy 


VLn(y)dy 


f"  1  d¥im(y),.  .  .  f"  1 
Jo  e(Y)  dv  Jo  e(y) 


1  d\j/4n,(y) 


‘-"i' 


e(y)  d  y 

1  d\|/im(y) 

e(y)  d  y 


e(y)  d  y 


¥4n(y)dy 


't'i„(y)dy  =1 


1  dy4ni(y) 

0  E(y)  d  y 


Vin(y)dy 


.  r^3.  .  P^3. 

i2mn=  vim(y)¥in(y)dy  =  1  V4in(y)¥in(y  M  y 

JO  •'0 

.hi  f‘’‘^  1  dVimCy)-  ^ 

•Mo 

h  f^^- 

l4mn=  ¥im(y)¥in(y)dy 
Jo 

f  ^  ~  -  f  *•  - 

lLn=  ¥lm(y)¥ln(y)dy  =J  ¥4m(y)¥4n(y)dy 
JO  *0 


^6mn' 


f"j_d^,„(y)-  r-L 

Jo  e(vl  dv  ^  Joey) 


Jo  e(y)  dy 


1  d\|/4n,(y)- 

■— — - ¥4n(y)dy 

e(y)  d  y 
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where  i=2,3,  for  regions  2  and  3,  correspondingly.  Thus,  an 

identical  set  of  equations  as  those  derived  for  h<y<h3  may  be 

derived  for  the  case  h3<y<h.  For  this  case,  it  is  sufficient  to 
substitute  i=3  instead  of  2  in  the  superscript  in  the 

orthogonalization  integrals,  and  substitute  3  for  2  in  the  subscript 
of  the  remaining  functions. 

By  developing  the  above  equations  further,  the  final 

form  of  the  set  of  eigenvalue  equations  is  derived  for  0<ySh3: 


Ml 


Ml  Th2  -e2 
^  ^3pTr  Imp 

^  2  2  e2 

(kz+k^,p)lfpp 


cos  y|  \  V  A  I 

sinf(klxmWl)  >Ain,+ 


M, 


V  1*^2 

Lj  ^2mri 


m=l 


S  2 

C  ]s°n<^lxmWl)+^"-%'^gQg(kixmWl)k2xTfOt[k2xn(w  >A„ 


(k^kLm)xh2 

Zj  .  -1  ^20111.^7- 


m=l 


2. 7-2  .  “““sinf[k2xn(W2-Wi)] 


A4m=0 


(3.1.19.a) 


(kz+k2xn) 


,  (k|*L^COS  COSlk2.,^W;-Wi)l 

1  jcos'*‘™  sinr*'’'"'  ''sinnk2,„(W2-w,)]  (  >' 


rono. 


'H  l3L!.!!f(klxmWl)-(kz+kLm)| 


m=l 


‘COS 


M,  Te2  .h2 
^  •*-4pn^2mp 

P-'  (k"+kLp)l5pp 


sinf  V  I  » 

QQg^kjxm'^l)  ^Ainf 


y‘  1  ^ 

(kz+kLn)  inf[k2xn(WrW7)]  "" 


(3.1.19.b) 
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m=l  (k^+kixn) 


-e2  COS,,  .. 

1 1  1  xm^  l)^  1  in+ 


if-  I 


M, 

^  ^Imn 
m=l 


\  (k^+kLn> 


COs[k2xn(W2-Wi)]+jk4^^sinf[k2xn(W2-Wi)]  )A4^- 


\, 


M, 


i  V  .-  -1|L  >sinflk2„(wrw  i)l  A4„=0 


M, 


e2  ,h2 


n;j2 


coM-o 


'm=l 


p-‘  (kf+kLpjijpp  j 


Asii 


M, 

s 


2  ~2 

(k^+kij^ni) 

(k^kLn) 


.h2  sinf  r 
^"’"cos  ^ 


Ixm 


iWi)Ai 


m 


(3.1.19.C) 


iSpjfn^p 

(kz+kLp)f2pp 


****  h2  ^ 

sinf[k2xn('^2-'^l)]-I  lmnSinfIk2xn(w2-w  i)l 


h2 

2mr' 


m=l 


(^^z+t^4xm) 

(k^kL,^ 


COS[k2xn(W2-W  i)]  +jk4xn^inf[k2xn(W2-W  j)] 


(3.1.19.d) 

As  before,  substituting  the  corresponding  variables  for 
region  3  will  yield  the  corresponding  equations  for  h3<y<h. 

The  above  set  of  equations  may  be  expressed  as  a 
homogeneous  matrix  equation  [A]x=0,  where  [A]  is  the  matrix  of 

the  coefficients  of  Ai^.  Ajm’  A4jn»  and  A4jn  x  is  the  vector 

containing  Aim,  Ai^,  A4m,  andA4in.  This  homogeneous  equation 
has  a  nontrivial  solution  only  for  the  case  det(A  )=0.  The 
coefficients  of  the  matrix  [A]  have  only  one  unknown  variable,  k^. 
so  solving  the  equation  det(A  )=0  yields  the  solution  for  the 
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propagation  constant.  By  substituting  this  value  of  back  into 
the  set  of  equations  [A]x=0,  an  underdetermined  system  is 

derived  for  determining  the  coefficients  Ajm.  Aijj,,andA4ni  of 

the  electric  and  magnetic  fields  of  the  structure. 


After  solving  for  the  coefficients  of  the  fields  in 
regions  1  and  4,  the  remaining  electric  and  magnetic  field 
coefficients  are  determined  by: 


^  ■  ^iym)  ^  2mm 

kj  V  sinf  ~  .T 

- L  l4pm^^  (klxpWi)Aip+ 

«liop=l 


Ip’ 


k,  ^>..2.-2  sinf~2 


■5:(k^klyp)*'"'(kupw,);  I 

n-1  hi  2-2 

'  l4.l(k  -kiyO 


(3.1.20.a) 


(3. 1.20. b) 


(3.1.20.C) 
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hi 


4mm 

(Op=i  Sini 


1 


,  ,  ei  2  2 

'='l2ii(k  -kiyO 


>Aip+ 


Ml  hi 

2i5;. 

p=i 


1 

■k  1 X 


py 


cos 

sin 


(kixpWi)Aip 


(3.1.20.d) 


The  coefficients  determined  through  the  matrix 
equation  are  determined  within  a  multiplicative  constant,  so  in 
fact  the  electric  and  magnetic  field  distribution  in  each  region,  if 
not  their  true  amplitude,  may  be  derived. 


3.2.  Characteristic  Impedance 


Theoretically,  three  different  expressions  may  be  used 
to  calculate  the  characteristic  impedance  of  any  device.  When 
applied  to  the  coupled  microstrip  line,  these  three  expressions 
define  the  relationship  between  the  characteristic  impedance  of 
the  coupled  microstrip  line,  the  time-averaged  power  flow  in  the 
strip,  the  complex  voltage  of  the  strip  calculated  at  the  center  of 
the  strip,  and  the  current  of  the  strip. 

The  power-voltage  relationship  gives  the  characteristic 
impedance  of  the  coupled  microstrip  line  as  a  function  of  the 
time-averaged  power  flow  in  the  strip,  Pavg»  and  the  complex 
voltage  V  of  the  strip,  calculated  at  the  center  of  the  strip.  In  the 
case  of  the  coupled  microstrip  line,  the  characteristic  impedance 
would  then  be: 


Z 


0  = 


(3.2.1) 


The  power-current  relationship  gives  the 
characteristic  impedance  of  the  coupled  microstrip  line  as  a 
function  of  Pavg»  and  the  complex  current  I  on  the  strip,  and  is 
given  by: 


(3.2.2) 


The  voltage-current  relationship  gives  the 
characteristic  impedance  of  the  coupled  microstrip  line  as  a 
function  of  the  complex  voltage,  V,  and  the  complex  current,  I,  of 
the  strip,  and  is  given  by: 


Although  experiments  have  shown  the  power-current 
calculation  to  be  more  valid  in  the  case  of  the  microstrip  line,  all 
three  calculations  will  be  presented  in  this  work. 

The  power,  P^vg  is  calculated  using  the  Poynting  vector, 
which  is  defined  as; 

P=YReJ  J  {ExH*aj}dxdy 

s 

Using  the  rectangular  compoents  of  the  electric  and  magnetic 
fields,  the  Poynting  vector  may  be  written  as; 

P=iRejJ(E,H;-E,H;)dxdy 

S 

Using  the  relationships  in  (3.1.8)  through  (3.1.11)  to 
express  the  fields  in  terms  of  the  scalar  potential  functions  in  each 
region,  a  Poynting  vector  may  be  calculated  for  each  region  of  the 
coupled  microstrip  line. 

The  power  in  region  1  is  thus  calculated  as; 
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.h  .w 


Pi=~ReJ^dyJ^  dx  (EiyH*x-EixH*y) 


^1  _V  T« 

1  2  ~9  h  /  2  _  _ 

=~7-  2w  2^  "‘^lyn^^SmnC  ,  (0,'VlMiy,m^O,k^^^,0)A^^Air,+ 


^  M^O  m=l  n=l 


1  I 


k  u  K  _  *  2 

(Op.Q  ~^lym)^5min  (0*Wi,kjxni»0,kixni»0)A  + 

m=l  1 

k 

^  (k^“k^ynj)l5jj,jri  (0,Wi,k2xni.0»kjxm»0)A.?n,  + 


Ml 


1  I 


,  ^lyn)^6nm  (  _  ) 


^  M-0  m=l  n=l 


-k2  U 
*^lxn  2 


where  the  functions  I1-I3  are  defined  as: 

ul 

l2(ll,ul,a,b,c,d)=  cos[a(x-b)]cos[c(x-d)]dx 
Jii 

ul 

l2(ll,ul,a,b,c,d)=  sinf[a(x-b)]sinf[c(x-d)]dx 
Ji  I 

f“' 

13(11, ul,a,b,c,d)=  cosf[a(x-b)]sinf[c(x-d)]dx 
Ji  1 

and  the  other  parameters  have  been  defined  in  section  3.1, 
The  power  in  regions  2  and  3  are  calculated  as: 
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1  /-h  [‘'^2 

Mj  Mj 

=-4-  SX<'‘  ^-^iy» ) '  3mn^MT  <mn- 
«  ^^0m=ln=l 

k 

m=l 

k  ^ 

-ji  X  (''^-'‘?yn,)l!i.m^NT3i,„+ 

m=l 

Mj  Mj 

4XS<>=^-‘'?ym)l4L.lNT4L 

m=  1 n= 1 


where  i=2  or  3,  for  region  2  or  3,  correspondingly.  The  functions 
INTI  through  XNT4  are  defined  as: 

T-NTlj^n^^^im^ixm^  1  ’  "^2’ ^ixn>  1  1 ) '*' 

®  im^  in^  1  1 ’''^2»^ixm»''^  1  ’  ^ixn’'^  1 )  ~ 

^im^ixm®  inl2(Wi,W2,ki^jn,Wi,kixn,Wi)  + 
BiniBinl3(Wi,W2,ki;,niyWi,kixn,Wj) 

TNT.^n~  ■^ini'^in^  1  1  ’"^2’  ^ixm>'^  1  >  ^ixn*'^  1  )'^ 

AimBinl3(Wi,W2,kixm>Wi,kixn.Wi)  + 

Binr^ in^3(''^  l»''^2’^ixn»''^l’^ixm’'^l)'*' 


®  im^  1  »''^2»^ixn’''^  1 ’^ixm*'^  l)"^ 

^im®  l’'^2»^ixm»'^l»^ixn*''^ 

BimBinl2(Wl,W2,kixin.Wi,kixn,Wi) 

XNT4njjj=— AjjjjAjjjk  jjjjjl3(W  J  *W2,kjxm»W  i*kjxn»Wj)  + 

^imBin^lC'^  l»'^2’^ixm’'^l  »^ixn»'^  l)“ 
Bim^in^ixn^2('^  l’''^2»^ixm»'^  1  »^ixn»'^l)'*‘ 
BiitiBin^S^'^  l»''^2’^ixn»''^l»^ixm»''^l) 

The  power  in  region  4  is: 
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The  total  power  in  the  device  is  equal  to  the  sum  of 
the  powers  calculated  for  each  region  of  the  line.  It  still  remains, 
however,  for  the  voltage  across  the  strip  to  be  calculated.  In  this 
case,  the  assumption  is  made  that  this  voltage  may  be 


approximated  by  the  voltage  at  one  point,  the  center  point,  of  the 
strip.  Thus,  the  voltage  may  be  defined  as: 


where: 

E2ix=2!ap2.yj=_J^£  (k^-k2ym)¥2m(y)- 

|A2mCOs[k2xm(w/2)]+B2mSinf[k2xni(w/2)]| 

Integrating  from  0  to  h3  yields: 

M2 

V=  ~^(k^-k2yni){A2mCOS[k2xm(w/2)]+B2mSinf[k2xni(w/2)]}- 
ni=i 

A  B 

-  — sinf(k2ymihi)+~sinf[k2ym2(h2-hi)]+ 


C  D 

- ; - {  1  -COS[k2ym2(h2-hl)] }+— COS[k2yn73(h3-h2)] 

^2^2ym2  J 

where  A,  B,  C,  and  D  are  the  expansion  coefficients  defined  in  eqs. 
(2.2.6.a-2,2.6.d)  for  the  TMy  potential  for  the  three-layer  parallel- 
plate  waveguide. 

The  current,  I,  on  the  strip  is  calculated  from  the 
transverse  magnetic  field  on  the  strip.  For  a  strip  of  zero 
thickness,  the  current  I  is  given  by: 
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[=J  H 

•'Strip 


transverse 


ds= 


r'^2 

H2y(y=h3)dx  +  H3y(y=h3)dx  = 

•’Wj  •'Wj 


M, 


^{j-k2-d-[(A2m-sinf(k2xm-w)+2-B2n,-sinf^(k2xm-w/2)]+ 


m=l 


1 


■d-[A2ni-(l-cos(k2x-w))-B2m-sinf(k2x-w)}- 


M3 

^{j-k^-cos[k3yni-(h-h3)]-[A3n,-sinf(k3xni-w)+ 

m=l 


2-B3n,-sinf^(k3xn,'W/2)+-r^ - C0S[k3  n,-(h-h3)]- 

[A3n,‘(l-COS(k3xni-w))-B3m-sinf(k3xni-w)]} 

The  equations  are  now  all  in  place  to  complete  the 
analysis  of  the  three-layered  coupled  microstrip  line  through  the 
calculation  of  the  propagation  constant  and  the  characteristic 
impedance  of  the  line. 
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Chapter  4.  Numerical  Results 

The  formulas  derived  for  the  calculation  of  the 
propagation  constant  and  the  characteristic  impedance  of  the 
coupled  microstrip  line  were  implemented  in  Fortran  77.  The 
Fortran  program  was  run  in  double  precision  on  the  VAX  or  in 
single  precision  on  the  CRAY,  with  the  same  results.  Some 
characteristic  results  are  presented  in  this  chapter. 

4.1.  Convergence  Criteria 

The  exact  electric  and  magnetic  field  representation 
for  each  region  of  the  waveguide  described  in  this  thesis  was 
given  by  an  infinite  sum  of  expansion  functions,  which  were 
characteristic  for  each  region  of  the  device.  Due  to  subsequent 

truncation  of  each  infinite  sum,  the  propagation  constant  which 

was  derived  for  this  structure  was  not  exact,  and  the  error 

present  in  the  calculation  not  exactly  known.  However,  as  was 
presented  in  chapter  2,  care  was  taken  so  that  the  relative 
number  of  expansion  functions  retained  for  each  region  was  such 
that  the  truncation  of  the  infinite  series  would  lead  to  a 
convergent  system.  The  accuracy  of  the  results  was  determined 
via  convergen^.,e  plots,  such  as  those  shown  in  figs.  (4.1.1)  and 
(4.1.2),  In  these  plots,  the  normalized  propagation  constant  was 
plotted  against  the  number  of  modes,  M2,  retained  in  region  2  for 
the  calculation  of  that  particular  value  for  the  propagation 
constant.  As  M2  increases,  the  value  of  the  normalized 
propagation  constant  converges  to  a  particular  value.  So,  after 
setting  a  certain  error  bracket  as  acceptable,  the  minimum  M2 
required  to  calculate  the  propagation  constant  within  that  error 
bracket  is  found  from  the  convergence  plots,  and  the 

corresponding  normalized  propagation  constant  is  accepted  as  the 
solution. 
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This  test  cannot,  however,  be  conducted  for  each  data 
point,  as  the  time  taken  to  do  this  would  be  prohibitive.  The 
procedure  followed  is  to  calculate  the  smallest  Mj  which  satisfies 
the  convergence  conditions  at  the  extreme  data  points.  These 
points  in  this  particular  case  would  be  smalles  and  largest 
structure  size,  and  lowest  and  highest  frequency.  It  is  safe  to 
accept  that  this  M2  will  then  satisfy  the  same  convergence  criteria 
for  any  other  point  within  the  boundaries  set  by  these  extreme 
points. 

The  error  bracket  for  the  convergence  plots  are  0.5% 
for  the  propagation  constant,  and  1.5%  for  the  characteristic 
impedance. 


Figure  4.1.1.  Convergence  plot  of  the  even  mode 
propagation  constant  versus  the  number  of  modes  retained  in 
region  2. 


Another  important  numerical  consideration,  which  has 
been  touched  only  briefly  so  far,  is  the  height  at  which  the  top 
plate  must  be  placed,  such  that  it  does  not  perturb  the  solution.  It 


56 


was  suggested  in  [19]  that  placing  the  top  plate  at  a  height  double 
the  height  at  which  the  metal  strips  are  placed,  would  be  enough 
to  ensure  an  unperturbed  solution.  This  is  the  height  used  for  all 
the  calculations.  The  same  program  was  run  for  greater  heights, 
but  it  was  deduced  that  this  only  added  unnecessarily  to  the 
computation  time,  without  producing  more  accurate  results. 


Figure  4.1.2.  Convergence  plot  of  the  odd  mode 
propagation  constant  versus  the  number  of  modes  retained  in 
region  2. 


4.2.  Program  Verification 

The  accuracy  of  the  programs  was  verified  by 
selecting  the  permittivities  of  the  dielectrics  and  the  geometry  of 
the  structure  such  that  it  would  be  reduced  to  structures  for 
which  results  were  available.  Both  the  propagation  constant 
calculation  and  the  characteristic  impedance  calculations  for  the 
single-layer,  coupled  microstrip  line  were  verified  against  the 
results  obtained  by  [2].  The  calculation  of  the  propagation 
constant  for  the  single-line  case  was  achieved  by  placing  the 
striplines  at  a  sufficiently  large  distance  from  each  other.  The 
results  obtained  for  this  case  for  the  single-layer  case  were 
verified  against  [2],  and  for  the  multiple-layer  case  against  [19  (p. 
57)].  The  calculation  of  the  propagation  constant  for  the  three- 
layer  coupled  microstrip  line  was  verified  with  [19  (p.  57)]. 

A  comparison  of  the  mode-matching  technique  versus 
the  spectral-domain  technique  was  achieved  through  the 
cooperation  of  Mr,Yu-De  Lin,  whose  analysis  may  be  found  in  [18]. 
The  two  different  programs  were  run  for  the  same  coupled 
microstrip  line,  and  the  results  are  presented  in  the  graphs  in  fig. 
4.2.1.  The  results  obtained  from  the  mode  matching  method  agree 
quite  well  with  those  obtained  through  application  of  the  spectral- 
domain  technique. 
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Figure  4.2.1.  Propagation  constant  and  characteristic 
impedance  (power-current  relationship)  verification  for  the 
coupled  microstrip  line. 


4.3.  Results-Design  Charts 


The  program  was  run  for  several  cases  of  coupled 
microstrip  lines.  Each  graph  shows  a  family  of  curves,  calculated 
for  a  specific  microstrip  width  and  dielectric  layering,  but  for 
different  strip  separations,  or  for  a  variable  dielectric  constant  in 
the  conducting  layer.  In  each  case,  for  wider  strip  separation,  the 
results  converge  to  the  case  of  a  single  microstrip  line.  Also,  in 
the  case  of  a  variable  dielectric  constant  value  in  the  conducting 
layer,  the  results  converge  to  the  single-layer  solution  when  the 
dielectric  constant  of  the  conducting  layer  approaches  the  value  of 
the  dielectric  constant  of  the  insulating  layers. 

Many  curves  were  graphed  versus  a  normalized 
frequency,  where  the  normalization  took  the  form  27th  3 Aq,  where  Xq 

is  the  free-space  wavelength.  The  characteristic  impedances, 
however, which  are  calculated  using  the  power-voltage,  Zo=V.V*/P, 
and  the  power-current,  Zo=P/I.I*,  expressions,  exhibit  discontinuities 
at  regular  intervals,  which  seem  to  correspond  approximately  to  the 
frequency  where  h3=nX(j/4,  where  n=l,2,3,etc.,  and  Xj  is  the  free- 
space  wavelength  in  a  medium  of  relative  dielectric  constant  equal  to 
the  relative  dielectric  constants  of  the  layered  medium.  This 
behavior  is  more  evident  where  the  characteristic  impedance  is 
plotted  against  h3/Xjj.  As  the  relative  dielectric  constant  of  the 
middle  layer  is  increased  with  respect  to  the  outer  layers,  this 
discontinuity  becomes  more  pronounced,  as  may  be  witnessed  in  figs. 
4.3.2(b),  and  (c).  This  behavior  seems  to  signify  the  "turning  on"  of 
higher-order  modes,  which  introduce  a  numerical  instability  around 
that  region  of  frequencies.  In  figs.  4.3.6(b),  and  (c),  the  characteristic 
impedances  are  plotted  over  a  wider  range  of  normalized 
frequencies,  and  the  almost  "periodic"  behavior  of  these 
discontinuities  is  more  evident.  The  calculation  of  Zq  using  Zo=V/I 

does  not  exhibit  any  discontinuous  behavior. 
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Normalized  Frequency 


Figure  4.3.4.  Zo=V/I  versus  27ih3Ao. 
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Figure  4.3.4rb'/.  Zo=P/n*  versus  h3Ad. 
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Figure  4,3.4(c).  Zo=VV*/P  versus  h3Ajj- 


Normalized  Frequency  (h  3  /  X,  d  ) 


Figure  4.3.4(d).  Zo=V/I  versus  h3Atj. 
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4.3.6.  Observation  of  the  behavior  of  the 
discontinuities  in  the  characteristic  impedance 
calculation  over  a  wide  range  of  frequencies. 


Normalized  Frequency  (2k  h  3  /  A.  q  ) 


Figure  4.3.6(a).  kz/kn  versus  27th3/Xo. 
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Figure  4.3,6(c).  Zo=VV*/P  versus 
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Figure  4.3.6(d).  Zo=V/I  versus  27rh3Ao- 


Chapter  5.  Conclusions 

In  this  thesis,  it  was  shown  how  a  quasi-planar 
structure  like  the  coupled  microstrip  line  on  a  layered  dielectric 

substrate,  could  be  analyzed  using  a  variation  of  the  mode 
matching  technique. 

Mode  matching  is  a  powerful  tool  for  analyzing  planar 
and  even  quasi-planar  structures.  The  results  obtained  through 
mode  matching  have  a  high  degree  of  accuracy,  especially  the 
propagation  constant.  This  technique,  however,  is  unsuitable  for 
CAD  packages,  as  it  is  very  inefficient.  It  is  also  susceptible  to 
numerical  errors,  as  the  system  which  is  solved  to  calculate  the 
field  coefficients  is  defined  by  a  singular  matrix.  It  should  be 
noted,  as  a  matter  of  fact,  that  this  singularity  condition  has 
actually  been  imposed  on  the  matrix,  in  order  to  solve  for  the 

propagation  constant. 

The  behavior  of  the  even  and  odd  mode  propagation 
constants  of  the  coupled  microstrip  line  on  a  layered  substrate 

does  not  differ  from  that  of  its  single  substrate  counterpart, 
except  for  a  slight  increase  of  k^/ko  with  an  increase  in  the 

conducting  layer  dielectric  constant  with  respect  to  the  insulating 
layer  dielectric  constants,  as  can  be  seen  in  fig.  4.3.2(a).  As  far  as 
the  impedance  is  concerned,  it  does  not  seem  to  change  noticeably 
with  the  change  in  dielectric  constants.  The  change  becomes  more 
noticeable,  however,  at  higher  frequencies.  The  effect  of  the 
discontinuity  becomes  much  more  noticeable  for  a  higher 
conducting  layer  dielectric  constant.  Its  position  does  not  change 
much,  however,  as  the  values  used  for  the  dielectric  constants 

were  of  the  same  order  of  magnitude. 

The  impedance  curves  seem  to  be  relatively  flat 
between  singularities.  The  discontinuities  themselves  follow  a 
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very  consistent  pattern  which  seems  to  be  dictated  by  the 
relationship  used  to  calculate  the  impedance  (power-voltage,  or 
power-current).  Irrespective  of  the  relationship  used  to  calculate 
Zq,  however,  the  first  singularity  occurs  at  the  same  normalized 
frequency,  h3/>.d=0.25.  This  seems  to  be  indicative  of  a  cutoff 
frequency,  above  which  higher  order  modes  start  propagating, 
and  around  which  the  calculations  of  the  characteristic  impedance 
become  numerically  unstable.  It  is  not  practical,  in  any  case,  to 
fabricate  microstrip  devices  of  the  order  of  magnitude  of  Xjj/4,  so 
only  the  portion  of  the  curve  before  the  discontinuity  occurs  is 
significant. 


Appendix  A.  Notation 

A  consistent  notation  scheme  is  used  throughout  the 
text  where  subscripts  and  superscripts  are  used  to  denote  region, 
direction,  dielectric,  mode  number,  mode  type,  or  to  differentiate 
between  two  similar  functions. 

A  variable  “z"  may  take  the  form  Zabm  where  "a" 
would  be  the  cartesian  coordinate  (x,  y,  or  z),  "b"  would  be  a 
region  number  (1,  2,  3,  or  4),  and  "m”  would  be  the  mode  number. 

A  function  "J"  may  take  the  form  where  "a" 

would  be  a  function  number,  "m",  and  "n"  would  be  mode 
numbers,  and  "c"  would  be  the  region  number;  c=2  means  the 
inner  product  was  performed  with  respect  to  a  normal  mode  of 
region  2,  and  c=3  means  the  inner  product  was  performed  with 
respect  to  a  normal  mode  of  region  3.  Superscript  "b"  is  used  only 
in  the  orthogonalization  integrals:  b=e  means  the  inner  product 
was  performed  with  respect  to  a  TM  mode,  and  b=h  means  the 
inner  product  was  performed  with  respect  to  a  TE  mode. 

The  following  general  notation  was  used  to  build  the 
desired  variables: 


t 


A  tilde  over  any  of  the  above  variables  would  denote 
TE  quantities,  whereas  no  tilde  would  denote  TM  quantities. 


Appendix  B.  Derivation  of  Field  Equations 


A  wave  is  a  field  that  is  a  function  of  both  time  and 
space.  Electric  and  magnetic  fields  that  vary  in  time  and  space 
are  governed  by  physical  laws  which  are  expressed  in  four 
equations,  known  as  Maxwell's  equations.  For  a  wave  travelling 
in  a  medium  characterized  by  a  certain  permittivity,  e,  and  a 
certain  permeability,  p.  These  equations  are: 

a  H 


VxE  =  -ii 


at 


VxH  =  e^  +  J 

a  t 


V  .  D  =  p 

V  .  B  =  0 


(B.l) 


The  electromagnetic  field  equations  above  are 
expressed  in  terms  of  six  quantities; 

E,  the  electric  intensity  (in  volts  per  meter) 

H,  the  magnetic  intensity  (in  amperes  per  meter) 

D,  the  electric  flux  density  (in  coulombs  per  meter) 

B,  the  magnetic  flux  density  (in  coulombs  per  square  meter) 

J,  the  electric  current  density  (in  amperes  per  square  meter) 

p,  the  electric  charge  density  (in  coulombs  per  square  meter) 

The  boldface  script  is  used  to  denote  complex  quantities. 

The  ultimate  sources  of  an  electromagnetic  field  are 
the  current  J,  and  the  charge  p. 


The  continuity  equation,  which  is  based  on  the 
principle  of  conservation  of  charge,  is  implicit  in  equations  (B.l), 
and  simply  states  that; 


J=-i£ 

at 


(B.2) 


Maxwell's  equations  are  complemented  by  the  so- 
called  constitutive  relationships,  which  incorporate  the 
characteristics  of  the  medium  in  which  the  field  exists.  These 
equations  define  the  electric  flux  density,  D,  the  magnetic  flux 
density,  B,  and  the  current  density,  J,  with  respect  to  the  electric, 
E,  and  magnetic,  H,  intensities. 


D=D(E,H) 

B=B(E,H) 

J=J(E,H) 


(B.3) 


Maxwell’s  equations  along  with  the  constitutive 
relationships  serve  to  fully  describe  a  wave  travelling  in  a  known 
medium. 


In  a  source  free,  linear  medium,  the  constitutive 
relationships  take  the  form: 

D  =  8E 
B  =  ^H 

(B.4) 

Here,  £  and  }l  are  constants,  where  e  is  the  capacitivity 
or  permittivity,  and  jj.  is  the  permeabillity  of  the  medium.  Using 
Eg  and  )lo  to  denote  the  corresponding  variables  in  vacuum,  for  a 
perfect  dielectric  (<7=0),  one  has  e=ereo»  where  is  the  dielectric 
constant,  or  the  relative  capacitivity,  of  the  medium,  and  li=j4o  for 
most  linear  matter. 
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The  present  analysis  will  be  concerned  only  with 
source  free,  linear  problems,  where  the  wave  has  a  steady-state 
sinusoidal  time  dependence.  In  this  case,  the  complex  field 
equations  read: 


-  VxE=z(co)H 
VxH  =  y(0L))E 


V.D  =  0 
V.B  =  0 


(B.5) 


A  A 

where  y(C0)=jC0£,  and  z(co)=ja)(io  in  nonmagnetic  material. 


The  above  representation  of  the  field  equations  gives 
rise  to  the  definition  of  the  parameter  k,  the  wavenumber  of  the 
medium.  The  wavenumber  is  defined  as: 


- ^ ^ - 

~y(co)  z({o) 


(B.6) 


The  physical  meaning  of  the  wavenumber  is  that  1/k 
is  the  velocity  of  propagation  of  an  electromagnetic  disturbance  in 
an  open  space  filled  with  perfect  dielectric  material  with 
permittivity  e  and  permeabillity  Po* 


Taking  the  curl  of  equations  (B.5),  and  using  the  above 
representation  of  k,  the  complex  wave  equations  become: 


VxVxE-k  E  =  0 
Vx  VxH-k^H  =  0 


(B.7) 


In  these  equations,  it  is  implicit  that: 

V.E  =  0  andV.H  =  0 


(B.8) 
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so  that  a  simplified  form  for  the  vector  wave  equations  may  be 
derived: 


and: 


V^E  +k^E  =0 


(B.9) 


V^H  +  k^H  =  0 


(B.IO) 


The  geometry  of  planar  structures  allows  us  to  work 
with  rectangular  cartesian  coordinates.  In  this  case,  the 
rectangular  components  of  E  and  H  satisfy  the  complex  scalar 
wave  equation: 


V^Ei  +  k^Ei=0 
V^H,+  k^Hi  =  0 


i  =  X  ,  y  ,  z 

(B.ll) 


To  construct  a  pliable  solution  to  the  above  equations, 
the  field  is  expressed  as  a  magnetic  vector  potential  A  and  an 
electric  vector  potential  F,  as  shown  below: 

E=-VxF+^VxVxA 

y 

E  =  VxA+^VxVxF 

^  (B.12) 


These  expressions  for  E  and  H  give  rise  to  a  very 
useful  classification  of  the  solutions  of  the  wave  equation.  In  this 
classification,  axial  uniformity  (the  cross-sectional  shapes  of  the 
waveguide  do  not  vary  in  the  direction  of  propagation)  is 
assumed.  In  addition,  this  classification  is  for  fields  conforming  to 
the  homogeneous  vector  Helmholtz  equations  (source  free 
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problem).  Propagation  is  assumed  to  be  in  the  z  direction,  and  the 
z  dependence  is  assumed  to  be  of  the  form  exp(+jpz).  The 
Helmholtz  equation  is  separable,  so  a  solution  of  scalar  Helmholtz 
equation  of  the  form  f(z)g(x,y)  is  sought. 


Under  this  classification,  the  vector  magnetic  and 
electric  potentials,  A  and  F  respectively,  may  be  assumed  to  be 
directed  along  one  coordinate  only. 


Choosing  A  =  UxV  where  xjr  is  a  scalar  wave  potential  (a 
solution  of  the  scalar  Helmholtz  equation),  will  yield  an 
electromagnetic  field  given  by: 


Ex 


Ez 


o 

II 

X 

X 

\d\  j 

_  1  d\ 

y  3x  dy 

>'  3z 

ro 

1 

It 

a=-^ 

y  dx  dz 

^  3y 

(B.13) 


This  choice  of  a  magnetic  vector  potential  will  yield 
Transverse  Electric  to  x  (TE*)  modes.  The  main  characteristic  of 
these  modes  is  that  Ez=0  and  Hz^tO.  All  TE^  field  components  may 
be  derived  from  the  axial  component  Hz  of  the  magnetic  field.  All 
field  components  may  be  derived  from  the  axial  component  Hz  of 
the  magnetic  field. 

Similarly,  choosing  F  =  u  x  V  will  yield  an 
electromagnetic  field  given  by  : 
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This  choice  of  an  electric  vector  potential  will  yield 
Transverse  Magnetic  to  x  (TM^ )  or  E  modes.  The  main 
characteristics  of  these  modes  are  that  : 

and  H2=0  (B,15) 

All  TM^  field  components  may  be  derived  from  the  axial 
component  Hz  of  the  magnetic  field. 


The  third  group  in  this  classification  are  the  TEM 
modes,  or  the  transverse  electromagnetic  waves.  The  waves 
belonging  to  this  classification  have  no  Ez  or  Hz  component.  In 

this  case,  the  electric  field  may  be  found  from  the  gradient  of  a 
scalar  function  \}/(x,y),  which  is  a  function  of  the  transverse 
components  only,  and  is  a  solution  of  the  two-dimensional  Laplace 
equation.  True  TEM  waves,  however,  will  be  found  to  occur  in 
very  few  cases  (i.e.,  free-space,  parallel -plate  waveguide,  etc.).  In 
addition,  true  TE  and  TM  modes  will  usually  not  be  sufficient  to 
satisfy  all  the  boundary  conditions  of  most  structures.  In  such 
cases,  however,  linear  combinations  of  TE  and  TM  modes  will 
provide  a  complete  and  general  solution. 

As  might  be  anticipated,  TE^  and  TM^  are  not  the  only 
existing  TE  and  TM  solutions  to  a  specific  problem.  By  assuming 
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the  magnetic  vector  potential  to  have  the  form  A=Uy\|/,  the 
following  equations  for  the  fields  wil  be  derived: 
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This  choice  of  A  will  yield  Transverse  Electric  to  y  (TEV)  modes, 
whose  main  characteristics  are  that  Ey^^O  and  Hy=0. 

Similarly,  taking  F=Uy  Vj/  will  yield: 
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This  choice  of  electric  vector  potential  will  yield  Transverse 
Magnetic  to  y  (TMy)  modes,  whose  main  characteristics  are  that 
Ey=0  and  Hy^^O. 

Similar  cases  may  be  derived  for  TE^  and  TM^. 
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